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We calculate analytically the spectra of plasma waves and electromagnetic waves (EMW) in 
metallic photonic crystal consisting of the parallel thin infinite metallic cylinders embedded in the 
dielectric media. The axes of metallic cylinders form a regular square lattice in a plane perpendicular 
to them. The metal inside the cylinders is assumed to be in the high frequency regime ujt S> 1, where 
r is the relaxation time. The proposed analytical theory is based upon small parameters / <C 1, 
where / is the volume fraction of the metal, and kR <C 1, where k is the wave vector and R is the 
radius of the cylinder. It is shown that there are five different branches of the EMW that cover all 
frequency range under consideration except one very small omnidirectional gap in the vicinity of the 
■ frequency of the surface plasmon. However, at some directions of propagation and polarizations the 

gap may be much larger. The reflection and refraction of the EMW is also considered. The general 
£N| . theory of refraction is proposed which is complicated by the spatial dispersion of the dielectric 

constant, and one particular geometry of the incident EMW is considered. 
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I. INTRODUCTION 



There has been a growing interest in recent years in the theoretical and experimental study of the photonic crystals, 
1-, 2- or 3-dimensional periodic structures that exhibit gaps for the propagation of electromagnetic waves (EMW) 
in different frequency ranges. Those gaps are similar to the Bragg's reflection of the X-rays The important 
• difference between the two effects is that in the X-ray theory one can use high frequency expansion of the dielectric 
constant to consider the interaction of the EMW with the media, while in the optical and microwave range this 
is impossible. The physical origin of the difficulties which may appear at low frequencies is the existence of other 
than EMW degrees of freedom. In dielectric photonic crystals the polaritons might be important Q while in the 
metallic photonic crystals (MPC) these degrees of freedom are provided by free electrons. There are different ways of 
classification of the MPC. The first criterion is dimensionality, they may be 1-dimensional 2-dimensional Q, and 
■ 3-dimensional ^ [|. The second criterion is the connectivity. The MPC may consist of the separated metallic islands 
fiof or it may include a connected metallic mesh . Both cases are considered theoretically in the review Jll|] . The 
third criterion reflects the frequency dispersion of the metallic conductivity. This dispersion is negledgible if ujt <C 1, 
where u> is the frequency of EMW and r is a relaxation time of electrons in the metal. In the opposite case, ujt 1, 
the dispersion of conductivity is very important, and the dielectric constant of the metal can be written in the form 

o : 2 

e p {u) = l-% (1) 



where uj 2 = Anne 2 /m is the plasma frequency of the metal, n - three-dimensional electron density, m is the effective 
mass of an electron. The low frequency case is considered in the majority of the works. Some works, however (see, 
say Ref. P,|l2|), devoted to the high frequency case. 

One of the most perspective methods of fabrication of the MPC is the infiltrating various metals and semi-metals 
into the opals using high pressure [p~3|Jl4| . In our numerical estimates below we mostly keep in mind these materials. 

In this work we consider a two dimensional MPC consisting of metallic circular cylinders in z direction, embedded 
into a dielectric media with the dielectric constant e m . The axes of the cylinders form a square lattice in the x-y 
plane. We assume a high frequency case so that dielectric constant of a metal is given by Eq. ([!]). 

The important applications of MPC are connected with the idea of creation a "transparent metal" , the material, 
which has a metallic type conductivity and at the same time is transparent at infrared or visible range. A regular 
bulk metal has the plasma frequency of the order of a few electron volts and it completely reflects light of the lower 
frequency. We show here that our system reminds a transparent metal even for the s-polarization (electric field is in z 
direction). It has a very narrow gap connected to a surface plasma wave and it is transparent at all other frequencies. 

A similar system has been considered theoretically by Kuzmiak et. al. Q. They assumed that the wave vector 
of EMW is perpendicular to the axes of the cylinders. Their computations show that for s-polarization the system 
completely reflects EMW with the frequency below U) V ^J //e m , where / is the volume fraction of the metal in the 
system. We confirm this result but we show that for the arbitrary angle of incidence the transparency is non zero at 
any low frequency. An important advantage of our method is that it is pure analytical. This gives us a deeper insight 
in the physics of the wave propagation, which is difficult for computational methods. Our method based upon small 
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parameters / < 1 and kR <§; 1, where R is the radius of the cylinders. For the case considered by Kuzmiak et. al. 
Pj our result for s-polarization can be obtained by a regular expansion of the system of equations derived in Ref. Q| 
with respect to parameters mentioned above. 

The main feature of our physical picture is a pretty complicated spectrum of the plasmons in the system, Slo (k) . Now 
we consider a few limiting cases to get a qualitative idea about this function. Gcnerically this spectrum comes from 



the one dimensional plasmon. The spectrum of the one dimensional plasmon for one cylinder only has a form 115 16 



" (fcz) = ^r l0g M?' (2) 

where R is the radius of the cylinder. This spectrum is almost linear and it starts from zero frequency. The electric 
field of the plasmon outside the cylinder decays exponentially as exp(— k z p), where p is the cylindrical coordinate. 

The physics is more sophisticated for a set of parallel cylinders which form a square lattice in the plane (x, y) with 
a lattice constant d. If the plasma wave is homogeneous in the (x, y) plane so that k x = k y = 0, and k z d <C 1 then 
electrons in a given cylinder feel electric field of many other neighboring cylinders. Then electrons in all cylinders 
vibrate with the same phase. Such a mode is equivalent to a three dimensional plasma wave with a frequency 
(jJ p \J f / e m - If k z d ^S> 1 there is no interaction between cylinders and one gets the one dimensional plasma wave. 
However its spectrum starts with the frequency to p ^J f / e m and there are no modes with a lower frequency. 

On the other hand, if k x — k y = Tr/d, the phases of ID plasmons in the neighboring cylinders are opposite, the 
electric fields created by all cylinders at a given one cancel each other. In this case the plasma waves in the system 
are still one dimensional and its spectrum starts with zero frequency. 

Plasma wave with a frequency f2o(k) can be excited by EMW which has a non zero E z component. This excitation 
is resonant if to of the EMW coincide with the corresponding value of plasma wave. Note, that the bulk plasma wave 
can not be excited by EMW because the first one is pure longitudinal while the latter is pure transverse. The plasma 
waves in our system are generated by electrons moving in the cylinders only, while the EMW propagates everywhere 
outside and inside the cylinders. Therefore in this case the EMW does excite plasma modes. 

There is another plasma mode in our system, which can be excited by electric field with non zero component 
perpendicular to the cylinders. This is a surface plasma mode with a frequency to s — to v j '\/e m + 1. This mode also 
plays the important role in the propagation of EMW. 

In this paper we show that there are five different branches of the EMW which cover all the frequency range with 
lot 3> 1 with one omnidirectional gap in the vicinity of the to s . Except of this gap there is another gap at lower 
frequencies which is not omnidirectional and which depends on the polarization of the EMW. The gap obtained by 
Kuzmiak et. al. Q is a particular case of such gap at k z = and s-polarization of EMW. 

Our paper is organized as follows. In section || we derive the spectrum of the plasma waves in the system. In 
section III we study the dielectric properties of the system. They are described by the tensor of dielectric constants. 
In the principal axes (x, y, z) two components of the tensor are equal to each other, namely e xx — e yy = e±_ and 
they depend on the frequency only. The longitudinal component e zz = ey is a function of to and k. The non-locality 
reflects the extra degrees of freedom and the long range Coulomb interaction. To calculate ey one should take into 



account deviation of acting field from average field. In section |V| we consider EMW in the system. They may be 
classified as ordinary and extraordinary waves by analogy with the optics of uniaxial dielectric crystals. However since 
the dielectric tensor is a function of to and k we get two branches for the ordinary wave and three branches for the 
extraordinary one. In section |v| reflection and refraction of EMW is considered. As in the case of optics of uniaxial 
crystals we obtain double refraction of EMW. 



II. SPECTRUM OF THE PLASMA WAVES 



The two-dimensional model considered in this paper consists of infinite metallic cylinders of a circular cross-section. 
The cylinders are parallel to the z axis and form a square lattice with the lattice constant d. Space between the 
cylinders is filled with a dielectric with the dielectric constant e m . 

To calculate plasma waves one should assume that the three-dimensional electron density in the I th cylinder with 
coordinates (xi, yi) has a modulation 

m = n + n 1 cos{k z z) e^^'+^V'-^ . (3) 
In this case electrostatic potential from the I th cylinder at a point with coordinates (x,y,z) has a form 
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*l{x,y,z) = AcosikzZy^+^-^Ko (k z y/(x - x{f + (y - Vl ) 2 ) , (4) 

where A = 2Trn'eR 2 /e m , Kq - modified Bessel function of 0-th order. 

The potential at the point (x, y, z) created by the whole system of cylinders can be written as 

*(x,y,z)= ®l( x > y> z) = Acos{k z z)e i{ - k * x+k ^-^y{x,y), (5) 

all wires 

where 



<p(x,y) = J2 e- l(K( - x - n - d)+k y (v ~ n y d]) K I k zy J(x ~ n x d) 2 + {y - n y d) 2 + R 2 J 

n x ,n v — — oo 



(6) 



Since R <C d the term R 2 under the square root is important for x — n x d = y — n y d = only. Note that <p(x, y) is a 
periodic function with the period d with respect to both arguments. 
Equation of motion for electrons in a cylinder is 

dV z 9$ 

where V z is a drift velocity, so that the current density j z = enV z . The continuity equation is 

Using Eqs. (||, [?], ||) one can obtain the plasma frequency of the system 

nl(k x ,k y ,k z ) = uj 2 ko f (k x , k y , k z ), (9) 

where 

„ fojld 2 „ / 
< = toT*» (10) 



Mk x ,k v ,k z )= J2' e« k - n ' +k » n >> d K (k z d^nl+nl) + K (k z R). (11) 

n x ,n y — — oo 

Here means that we should exclude the term n x = n y = from the summation. One can see that the function 
f2o is a periodic function of k x and k y with the periods equal to reciprocal vectors of the square lattice, but it is not 
periodic with respect to k z . 

In the case when k x — k y — and k z d <C lwe obtain that (/?o(0, 0, k z ) = 2ir/k 2 d 2 , thus the plasma frequency of the 
system f2g(0, 0, k z ) — u> 2 (f /e m ). In another limiting case k x = k y = and k z d ^> 1 the plasma frequency is given by 
Eq. (||) for the one dimensional plasmon. The same result is valid at any k z if k x = k y — ir/d. 

Figs. |l|(a) and 0(b) illustrate the behavior of the £lo(k x , k y ,k z ) as a function of k x ,k y ,k z . 

III. DIELECTRIC PROPERTIES OF THE SYSTEM 
A. Calculation of the e zz = ey component of the dielectric tensor 

Let us introduce an external electric field acting in z direction in a form 

£ t = £' ^{k x x+k y y+h z z—ut) 0-2) 

Now the equation of motion has a form 

m ~jT~ = e (£cxt + £)■ (13) 
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Here £ is an internal field in z direction given by 



<9$ ., 2vrei? 2 



£(x,y) =-jr- = -ik z ri e^ k " x+ ^ v+K '*- ut >(p(x ) y). (14) 

OZ Cm 



Using the continuity equation (0) and the results of Sec. || one gets 

, ienk. 



n m(^-n 2 (k x ,k y ,k z )) £ ' ext - (15) 
Now we can relate current density j z to the external field 

zenu 

3z k z m(^-n 2 (k x ,k y ,k z )) Cext - [0) 

To find effective macroscopic conductivity a one should relate average current density j z to the average electric 
field E by equation 

J z = <tE. (17) 

The field E can be found by the direct averaging 

E = £ cxt + -J2 [ [ £( x > v) dxdy, (18) 

« J Jceli 

while 

J z = fjz- (19) 
Using equations ( |l4]) and (|l5|) one can obtain a relation between the average and the external fields 



~ip(k x ,ky,k z ) = -jo / / ip(x,y) dxdy (21) 
" J JceU 



where 



After double summation and integration in Eq. (Ell) one can obtain 



Using Eqs. (E6[) , jl7|) , (fl9|) , (pp|) one can find cr. Since the dielectric constant 



where 



A-k/j 

e = e m + i (23) 

one can obtain for the e zz = en component of the dielectric tensor 

= e m - u2 _ n2{k P xM y (24) 

2 2 / 2 & 2 

^ (k x , k y , k z ) — Q a (k x , k y , k z ) oj p 1,2 T l2 T 1,2 ' (2^) 

e m K x ' K y ' K z 

Until now we assumed that relaxation time of the electrons r is infinite. To describe the decay of the EMW one 
should substitute lo 2 by ui(u> + ir^ 1 ) in Eq. (pi|). 

Dependence of en on the frequency at a fixed value of the wave vector k is shown at the Fig. || for the case k z ^ 0. 
One can see that at small frequency ey is positive. The dielectric constant e\\ becomes zero at uj = 0J\ = (fl 2 + fui 2 ) 1 / 2 . 
Note that the region of the negative en is not omnidirectional and it disappears at some values of k. 
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B. Calculation of the e xx — e yy — e± components of the dielectric tensor 

Assuming that / <C 1 one can apply the method for finding the dielectric constant of the granular mixture, described 
by Landau and Lifshitz (see Ref. jL7|, p. 45). One starts from the identity 

i J (X>(r) - e m £ (r))dV = D e m E, (26) 

where D and E are average induction and electric field respectively, V is the volume of the system. In the following 
calculations we assume that all fields are perpendicular to the cylinders. The dielectric constant is determined by 
relation D = ej_E. The integrand in the left side of Eq. (^6|) is non-zero only inside the cylinders. Since / is small 
one can assume that electric field acting on a cylinder is E. The field inside the cylinder Ej„ is constant and equal 
(see Ref. @, p. 43 ) 

E m = -^J^E, (27) 

tin * 

where internal dielectric constant ej n is given by Eq. ([!]). As the result the component of the dielectric tensor ej_ is 

(em — l)^ 2 + 

e±M = e m -2/ em m 2 §. (28) 

(e m + 1 )ui — LU p 

Note that the pole in Eq. (Bsh corresponds to a surface plasma wave propagating along the surface of the cylinder. 

This equation is written in the first approximation with respect to the small parameter /. This approximation 
can be improved by a two-dimensional version of Clausius-Mossotti procedure. One can show that effective field is 
E + 2-7rP, where P is the polarization. As a result one gets an extra term 2/w 2 in the denominator of Eq.(p8|) which 
can be neglected. 

Finally, we have found all diagonal components of the dielectric tensor e zz — eii and e xx = e yy = e±. Since z is 
chosen in the direction of the axes of the cylinders the dielectric tensor is diagonal at any choice of x and y. 

IV. SPECTRUM OF ELECTROMAGNETIC WAVES IN THE SYSTEM 

The Maxwell's equation for the EMW has a form (see Ref. Q, p. 316) 

n 2 E l - m(nkEk) = t ik E k , (29) 

where 

n = - k. (30) 

LU 

The spectrum of EMW can be found from the equation 

det\n 2 S ik - n t n k - e lk (n x ,n y ,n z )\ = 0. (31) 

Since the structure of the dielectric tensor is the same as for a uniaxial crystal one can get that Eq. (|3l]) falls off 
into two equations: 

n 2 = e x {uj) (32) 

and 

1 sin 2 8 cos 2 8 , . 

(33) 



n 2 e\\(n x ,n y ,n z ,u)) ej_(w) 
where 8 is the angle between wave vector k and z axis. Using analogy with the uniaxial dielectric crystals we call 



the solutions of Eq. (|32|) and Eq.(33) the ordinary and the extraordinary waves respectively. In the ordinary waves E z 



component of the electric field is zero. 
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In fact, our equations are much more difficult than in the regular optics of dielectric crystals because of the 
complicated dependences of the dielectric tensor on u and k. Due to these dependences there are two ordinary waves 
and three extraordinary waves. 

The frequencies of the two ordinary waves depend only on |fc| and can be obtained analytically 



e m (l + 2f)m 2 p + c 2 k 2 (e m + l) 
2e m (e m (l-2/) + l + 2/) 



J (e m (l + 2/)u; 2 + c 2 k 2 (e m + l)) 2 - Ac 2 k 2 cu 2 e m (e m (l - 2/) + 1 + 2/) 

± — (34) 

2e ro (e m (l - 2/) + 1 + 2/) V ; 

These solutions are shown in the Fig. ||. All the calculations shown in the plots below are made for / = 0.1 
and c/u>pd = 0.03. The lower branch of the ordinary wave tends to the frequency of the surface plasmon mode 
lo s = ujp/ V e m + 1 as ck/ujp tends to infinity. The upper branch starts with the frequency u> s (l + /e m )and transforms 
into a regular dispersion law for the EMW in the system with the dielectric constant e m (l — /). The spectrum of the 
ordinary wave has a gap of the width ft m uj s . As we show below, a part of this gap is covered by the extraordinary 
waves, however a small omnidirectional band gap still exist in this region. 

To find a dispersion relation for the extraordinary wave one should solve cubic equation for to 2 (Eq. J33|)). Since 
the analytical solution is not instructive we present only the results of numerical calculations. The component of the 
dielectric tensor cy depends on \k\,9, and the azimuthal angle ci>. 

Now we demonstrate the results of our calculations of the EMW spectra in the first Brillouin zone only. The EMW 
spectra in higher zones can be easily obtained but since we keep in mind small values of d they correspond to pretty 
high frequencies. The 3D plots Fig. ||(a) and (b) show three branches of extraordinary waves in coordinates ck/ui p 
and 9 from different points of view. The azimuthal angle <j> is chosen 45° for all plots below. Fig. || shows three 
extraordinary branches at 9 — 57°, 4> = 45°. Fig. || shows extraordinary branches at 9 = 0. One can see from Eq. (|3^ ) 
that at 9 = the extraordinary branches coincide with the ordinary ones. Fig. || shows extraordinary wave at 9 = ir/2. 
Note that there are singularities at 9 — and 9 = ir/2 (see Eq. j3^)), because at 9 = the zero of ey does not play 
any role. The same happens with the zero of e± at 9 = ir/2. Therefore Fig. || has only two branches. The lower 
extraordinary branch in Fig. || is the result of unification of the lower and the middle branches of the Fig. f|(a). They 
form one branch which coincide with the lower branch at small k and with the middle branch at large k. The same 
happens with the middle and upper branches in Fig. ^(b). They form one branch which coincide with the middle 
branch at small k and with the upper branch at large k. Since the lower branch is absent at 9 — ir/2, Fig. ^| has only 
one branch which coincides with the results of the Ref. 0. Note, that all singularities will be smeared out if one takes 
into account the finite value of r. 

The lower branch in Fig. [| starts with the linear dispersion law u = ck cos 9 j ^ft^ . At k z l/d this branch tends 
to the one dimensional plasmon with a dispersion law Eq. (||) , but its frequency remains higher than the frequency of 
the plasmon. 

The middle branch starts with the frequency u> p (f /em) 1 ' 2 and saturates at w s ( 1 + / sin 2 9(e m — 1)/ (e m + l)), which is 
above saturation value of the lower branch of the ordinary wave lu s . The upper branch starts at the w s (l + /e m ). Note, 
that the upper branch of the ordinary wave starts with the same frequency. It is easy to show that the omnidirectional 
gap is very narrow and equal to fiv s (e 2 n + l)/(e m + 1)- At all frequencies below and above this gap the system is 
transparent. The upper branch transforms into the regular photon with the dispersion cfc/e m (l — /). 



V. REFLECTION AND REFRACTION OF ELECTROMAGNETIC WAVE 

Fig. ^ shows complete spectrum of electromagnetic waves in the system, which includes both ordinary and extraor- 
dinary waves. One can see, that every frequency corresponds to two waves (ordinary and extraordinary) everywhere 
except the gap for the ordinary waves. Thus, in a general case, one should expect a double refraction similar to a 
regular optics of dielectric crystals. Our case is more difficult because of the ui and k dependence of the dielectric 
constant. 

We mention first that tangential component of the wave vector of incident, reflected and refracted waves should 
be the same. This follows from the very general statement (see Ref. [pi, p. 304) that the spatial variation and time 
variation of the electric and magnetic field vectors must be the same at the boundary. Suppose that the boundary is 
the plane £ = and the plane of the incidence is (£, 77). Then one gets 



G 



n„, = n lv = n 2r) = n vq = sin , (35) 

where symbols r, 1, 2, i label reflected, ordinary refracted, extraordinary refracted and incident wave respectively, Oi 
is the angle of incidence. 

The solution of the problem of reflection and refraction in the case of k dependent dielectric constant should consist 
of the following parts. The first part is to find and n 2 ^ which determine angles of refraction of ordinary and 
extraordinary waves 0\ and 2 respectively. For the ordinary wave refraction angle can be found using the SnelPs law 

sm#i = — — . (36) 



Here and below we assume that incident and reflected waves propagate in the vacuum. 

Due to spatial dispersion of the eii calculation of the refraction angle for extraordinary wave is more complicated. 
One should find a solution for the frequency of the proper branch of the extraordinary wave in a form u> = f(ii 2 )- 
This equation can be used to find n 2 £ at given and u. The refraction angle for the extraordinary wave can be 
found using relation tan 9 2 = n^jn^. 

The second part of the problem is to find electric field of reflected E r = (E rn , E r ^, E r ^), ordinary refracted Ei = 
(Ei v , E\q, E\^) and extraordinary refracted E2 = {E 2ri ,E 2 ^,E 2 ^) waves. For this purpose we should first find two 
relations between three components of vector E in each refracted wave from Eq. (^9|) . The other equations are the 
regular boundary condition on the surface including the continuity of the normal component of the displacement 
vector Di = e ik E k . 

Strictly speaking the last condition should be formulated in coordinate space rather than in k space Jig) ] . However 
it works in k space if kd 1. In general case one can use the continuity of the normal component of the Poynting 
vector, calculated far from the boundary. The scheme described above is valid when the boundary plane forms any 
angle with the axis of cylinders in our system. Note that if electric field in the incident wave is perpendicular to the 
cylinders, there is only ordinary refracted wave. In this case equations for all fields have a regular form (see Ref. p8| , 
p. 305). 

As an example we consider a case when the boundary plane is perpendicular to the cylinders (axis £ coincides with 
the z axis) and the plane of incidence includes z axis and the diagonal of the square lattice formed by the axes of 
the cylinders in the plane z = 0. Thus for all waves n x — n y — sin Oi /a/2- For the ordinary wave nf z — e± — sin 2 Oi. 
Suppose that we have found n 2z (u>) for the extraordinary wave. Note that in different frequency ranges n 2z {u>) is 
determined by different branches of extraordinary waves. Then for electric fields of the reflected and refracted waves 
one gets 

Erx = ^ _ _L_ j Eix + (^i - E iv , (37) 

E ry = ( - -J±-) E ix + ( 1 —) E iy , (38) 

\ A n lz + lj \ A n lz + lj 

( cos 2 6i tan 6 2 e± \ , 
E rz =[l-2 ±- - 4, 39 

Eix = ~Eiy — ——{Eix — Eiy), E\ z — 0, (40) 

7l\ z T" 1 

E 2x — E 2y — — 1 (E ix + E iy ), (41) 

cos 2 Oi tan 2 e± , , „ 

E 2z = 2 ^- i— ± -E lz , (42) 



sin6>i Ae\\(n X: ny,n 2z ) 



where 
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A = tan 60 sin t 



i e\\(n x ,ny,n 2z ) 

Reflection coefficient 1Z can be found using Eqs. d3HH) as 

R < r 



- e_L 



cos v t + n 2z 



K = 



Ei 



(43) 



(44) 



Fig. H shows the reflection coefficient in the case of normal incidence. For this case E z = in reflected and refracted 
waves and only the ordinary wave propagates through the system. To avoid singularities we assumed finite value of 
the relaxation time. One can see that the system is basically transparent except for the small gap in the spectrum of 
EMW in the vicinity of the surface plasmon. 



VI. CONCLUSION 



We have found that the two dimensional metallic mesh has two different plasma modes. One of them comes 
generically from one dimensional plasmon, which is produced by electrons moving along the thin metallic cylinder 
(wire) with both longitudinal and transverse electric fields outside the cylinder. It has a complicated dispersion law 
and its frequencies cover all the range from zero to the effective plasma frequency u> p y/ f /e m . Such a plasma mode 
can be excited by s-polarized EMW and it contributes to the longitudinal component of the dielectric tensor eii, which 
is a function of the wave vector k, because the plasma frequency depends on k. Another important plasma mode is 
a surface mode with the frequency u v j \Jt m + 1. This mode is almost independent of k. It is excited by p-polarized 
EMW and it contributes to the perpendicular component of the dielectric tensor ej_, which is fc-independent. 

These two plasma modes generate five different branches of the EMW, which cover all the frequency range except 
very narrow omnidirectional gap in the vicinity of the frequency of the surface plasmon. The modes can be classified 
as two ordinary waves (E z = 0) and three extraordinary waves [E z =/= 0). This classification comes from the uniaxial 
dielectric crystal, but the amount of modes in the metallic mesh is larger because the Fresnel equation becomes 
non-linear due to the frequency and the spatial dispersions of the dielectric tensor. 

The theory of reflection and refraction of the EMW in such a system becomes a special problem. We have described 
a general approach to this problem and solved it for the case of particular plane of incidence. The reflection coefficient 
has a maximum in the narrow frequency interval of the omnidirectional gap. 

Finally, this system is a good example of the "transparent metal" , which has a metallic conductivity in z-direction 
and is transparent almost in all frequency range where condition lot 3> 1 is fulfilled. 
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FIG. 1. Plasma frequency Qo(k) in the plane k x = (a) and in the plane k z d = 0.5 within one quarter of the first Brillouin 
zone (b) at the following values of parameters / = 0.1, c/u p d = 0.03. 
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FIG. 2. Dielectric constant en as a function of frequency at a fixed value of k with non zero k z . Relaxation time is taken to 
be infinite. 




FIG. 3. The spectra of the ordinary wave as given by Eq.(|34|) (solid lines). The lower dotted horizontal line is uj = lu s . The 
upper dotted line is u — w s (l + /f ra ). The dashed line shows the photon dispersion u — ck/e m (l — /). Note that the solid lines 
represent also the spectrum of the extraordinary waves at 6 = 0. / = 0.1, c/uj p d = 0.03 
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FIG. 4. Three branches of extraordinary waves in coordinates ck/u p and 9 from different points of view. / = 0.1, 
c/ujpd = 0.03. 
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FIG. 5. Three extraordinary branches at 8 — 57°, (f> = 45°. The lower dotted horizontal line is 
uj = u> s (l + /sin 2 #(e m — l)/(e m + 1)). The upper dotted line is uj = cu s (l + /em). The dashed line shows the photon 
dispersion uj — ck/e m (l — /). / = 0.1, c/uj p d = 0.03. 
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FIG. 6. The only extraordinary branch at 6 = 7r/2. The dashed line shows the photon dispersion uj — ck/e m (l — /). / = 0.1, 
c/uj p d = 0.03. 
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FIG. 7. Dispersion of ordinary (full line) and extraordinary wave (dashed line). This figure combines Fig. ^|and Fig. The 
dotted line shows the photon dispersion ui = cfc/e m (l — /). 
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FIG. 8. Reflection coefficient for the case of normal incidence for the following values of parameters lu p t = 10 (solid line), 
UpT = 5 (dotted line), e m = 1.5, / = 0.1. 
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